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Abstract 

Two  geometries  of  solid  oxide  fuel  cells  (SOFCs)  are  currently  under  development,  tubular  and  planar.  Both  types  of  cells  are  configured  in 
a  stack  using  an  interconnect,  which  electrically  connects  the  anode  of  one  cell  to  the  cathode  of  the  adjacent  cell,  while  also  physically 
isolating  fuel  from  oxidant.  Proper  design  of  the  interconnect  in  conjunction  with  single-cells  is  critical  to  minimizing  the  overall  stack 
resistance.  This  work  quantitatively  examines  the  dependence  of  total  SOFC  stack  resistance  as  a  function  of  interconnect  contact  spacing, 
interconnect  contact  area,  cathode  thickness,  electrolyte  thickness,  anode  thickness,  and  transport  properties  associated  with  each  region  and 
at  interfaces  (charge  transfer  resistance).  Both  one-dimensional  (channels)  and  two-dimensional  (dimples)  symmetries  of  interconnect 
geometry  are  analyzed  for  planar  cells.  Analytical  expressions  are  derived  for  the  area-specific  resistance  (ASR)  of  a  repeat  unit  consisting  of  a 
cell  and  interconnect,  for  both  geometries,  as  a  function  of  cell  parameters,  interconnect  contact  area  and  interconnect  contact  spacing.  It  is 
found  that  the  one-dimensional  interconnect  symmetry  leads  to  lower  values  of  stack-repeat  unit  area- specific  resistance  (ASR)  than  the  two- 
dimensional  symmetry.  Thus,  based  on  the  analysis  presented  here  the  one-dimensional  interconnect  geometry  is  preferred  over  the  two- 
dimensional  one. 

©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

The  performance  of  solid  oxide  fuel  cells  (SOFCs)  at  the 
single-cell  level  has  been  improved  dramatically  through  the 
use  of  thin  electrolytes  and  porous  composite  electrodes. 
With  techniques  such  as  dip-coating,  slurry  coating,  tape 
casting,  and  electrochemical  vapor  deposition  (EVD),  it  is 
possible  to  fabricate  gas-tight  yttria- stabilized  zirconia 
(YSZ)  electrolytes  of  thickness  in  the  range  of  5-30  pm, 
which  have  a  very  low  area- specific  resistance  at  tempera¬ 
tures  >700  °C  [1-3].  Porous  composite  electrodes,  which 
are  mixtures  of  electrocatalyst  and  electrolyte  materials,  or 
single-phase,  mixed  ionic-electronic  conducting  (MIEC) 
materials,  have  been  shown  to  significantly  reduce  the 
area-specific  resistance  associated  with  charge  transfer  (acti¬ 
vation)  reactions,  both  theoretically  and  experimentally 
[4-6] .  Theoretically,  the  effective  charge  transfer  resistance 
of  a  porous  composite  electrode  is  minimized  for 
relatively  thick  electrodes  (several  micrometers)  with  a  fine 
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micro  structure,  in  the  immediate  vicinity  of  the  electrolyte 
[3,6].  Single-cells  consisting  of:  (1)  a  porous  Ni  +  YSZ  anode 
of  thickness  ^500-1000  pm,  (2)  a  dense  electrolyte  of  YSZ 
of  thickness  ~5— 15  pm,  and  (3)  a  porous  Sr-doped  LaMn03 
(LSM)  +  YSZ  cathode  of  thickness  ~100  pm  typically 
exhibit  power  densities  (area-specific  resistances)  of 
~0.8  W/cm2  (-0.3 1  Q  cm2)  and  -1.5  W/cm2  (0.17  Q.  cm2) 
at  700  and  800  °C,  respectively.  Both  the  cathode  and  the 
anode  are  often  compositionally  and  microstructurally  graded 
through  the  thickness.  Near  the  electrolyte/electrode  inter¬ 
face,  the  electrode  microstructure  is  generally  fine,  and  this 
region  may  exhibit  mixed  ionic-electronic  conductivity. 
Away  from  the  interface,  the  electrode  usually  has  high 
porosity,  the  microstructure  is  usually  coarse,  and  the  con¬ 
duction  mechanism  is  predominantly  electronic. 

For  the  purpose  of  practical  power  generation,  several 
single-cells  must  be  connected  in  a  stack  to  achieve  a 
sufficiently  high  power,  and  build  the  desired  voltage.  In 
a  SOFC  stack,  the  single-cells  are  electrically  joined  with 
interconnects.  Ideally,  the  power  generating  capacity  of  the 
stack  should  be  the  sum  of  capacities  of  the  individual  cells. 
However,  there  are  resistances  associated  with  interconnect 
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Nomenclature 

E  Nernst  potential  (V) 

I\  total  current  for  the  one-dimensional 

interconnect  design  (A/cm2) 

7n  total  current  for  the  two-dimensional 

interconnect  design  (A/cm2) 

L  center-to-center  spacing  between  nearest 

interconnect  contacts  for  the  one-dimensional 
geometry  (cm) 

M  center-to-center  spacing  between  nearest  inter¬ 

connect  contacts  for  the  two-dimensional 
geometry  (cm) 

r  radial  position  measured  from  the  center  of  an 

interconnect  contact  (cm) 

r0  radius  of  the  interconnect  contact  with  the 

electrodes  (cm) 

R  area-corrected  spacing  between  interconnect 

contacts  for  the  two-dimensional  geometry 
when  circles  are  used  instead  of  hexagons  as 
the  in-plane  repeat  unit  (cm) 

Ri  overall  area- specific  resistance  of  a  repeat  unit 

for  the  one-dimensional  interconnect  design 
(cm2) 

Rn  overall  area- specific  resistance  of  a  repeat  unit 

for  the  two-dimensional  interconnect  design 
(Q  cm2) 

R*t  charge  transfer  resistance  at  the  interface  of 
the  electrolyte  and  a  composite  anode  (Q  cm2) 
Rlt(ty  intrinsic  charge  transfer  resistance  at  the 

interface  of  the  electrolyte  (e.g.  YSZ)  and 
the  electrocatalyst  (e.g.  Ni),  for  the  given 
particle  size  and  volume  fraction  (Q  cm2) 

Rcct  charge  transfer  resistance  at  the  interface  of 
the  electrolyte  and  a  composite  cathode 
(Q  cm2) 

Rcct(o^  intrinsic  charge  transfer  resistance  at  the 

interface  of  the  electrolyte  (e.g.  YSZ)  and 
the  electrocatalyst  (e.g.  LSM),  for  the  given 
particle  size  and  volume  fraction  (Q  cm2) 

4  thickness  of  the  anode  (cm) 

tc  thickness  of  the  cathode  (cm) 

4  thickness  of  the  electrolyte  (cm) 

vT(x)  electrostatic  potential  between  the  cathode 

and  the  anode  as  a  function  of  position  for 
the  one-dimensional  interconnect  design  for 
Vo  <  v  <  LI 2  (V) 

vn(r)  electrostatic  potential  between  the  cathode 

and  the  anode  as  a  function  of  position  for 
the  two-dimensional  interconnect  design  for 
r0  <  r  <  RI 2  (V) 

v  horizontal  position  measured  from  the  center 

of  an  interconnect  contact  (cm) 
v0  half-width  of  the  interconnect  contact  with  the 

electrodes  (cm) 


Greek  letters 

</>°  electrostatic  potential  maintained  at  the  cath¬ 
ode  interconnect  contact  relative  to  the  anode 
interconnect  contact  (V) 

4>i(x)  electrostatic  potential  as  a  function  of  position 
within  the  anode  for  the  one-dimensional 
interconnect  design  for  v  <  v  <  LI 2  (V) 

</)j  (x)  electrostatic  potential  as  a  function  of  position 
within  the  cathode  for  the  one-dimensional 
interconnect  design  for  x$  <  x  <  L/2  (V) 

</)jj  (r)  electrostatic  potential  as  a  function  of  position 
within  the  anode  for  the  two-dimensional 
interconnect  design  for  ro  <  r  <  RI2  (V) 
4>ii(r)  electrostatic  potential  as  a  function  of  position 
within  the  cathode  for  the  two-dimensional 
interconnect  design  for  ro  <  r  <  RI2  (V) 
pa  electronic  resistivity  of  the  anode  (Q  cm) 

pc  electronic  resistivity  of  the  cathode  (Q  cm) 

pe  ionic  resistivity  of  the  electrolyte  (Q  cm) 

Pg  ionic  resistivity  of  the  cathode  (e.g.  that  of 

YSZ  in  the  cathode)  (Q  cm) 


material  and  interconnect  design  parameters.  As  a  result,  the 
power  generating  capacity  of  a  stack  is  less  than  the  sum  of 
the  capacities  of  single-cells.  Therefore,  an  efficient  inter¬ 
connect  design  is  of  paramount  importance.  Some  of  the 
parameters  which  influence  the  resistance  contributed  by  the 
interconnect  can  be  deduced  by  considering  the  current  path 
for  electrons.  Fig.  1  shows  a  schematic  of  a  SOFC  cell- 
interconnect  repeat  unit  and  the  regions  of  contact  with 
interconnect.  Current  flows  from  the  electrocatalyst-elec¬ 
trolyte  interface,  through  the  electronically  conductive  elec¬ 
trocatalyst,  such  as  LSM  in  the  cathode,  to  interconnect. 
Usually,  the  resistance  of  interconnect  can  be  assumed  to  be 
small;  although  this  would  not  be  the  case  if  a  resistive  oxide 
scale  forms  on  the  interconnect.  The  overall  resistance  also 
depends  upon  the  spacing  between  the  interconnect  contacts 
due  to  the  resistances  associated  with  the  two  electrodes  and 
the  electrolyte  comprising  the  cell.  In  fact,  parameters 
related  to  the  electrodes,  which  are  not  formally  part  of 
interconnect,  influence  the  effective  stack  resistance.  For 
this  reason,  the  design  of  interconnect  is  influenced  by  the 
properties  of  the  single-cell.  That  is,  the  single-cell  and 
interconnect  should  be  designed  as  a  subsystem  and  not  as 
independent  entities. 

Several  numerical  models  for  solid  oxide  fuel  cell  stacks 
have  been  developed,  which  include  an  interconnect  [7]. 
Ahmed  et  al.  modeled  a  planar,  monolithic  cell  [8].  The 
model  addresses  the  effects  of  fuel  utilization  and  heat 
generation  on  the  average  current  density,  but  it  does  not 
include  resistance  to  distribution  of  current  in  or  over  the 
electrodes  parallel  to  the  electrolyte.  Bessette  et  al.  devel¬ 
oped  a  finite-element  model  of  the  tubular  Siemens-Wes- 
tinghouse  design,  which  treats  current  distribution  losses  in 
the  electrodes  with  a  resistor  network  [9].  Local  current 
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Fig.  1.  A  schematic  of  a  SOFC  in  series  contact  with  interconnect  of  one-dimensional  geometry. 


density  was  calculated  as  a  function  of  axial  position  and 
angle.  Local  current  density  was  noted  to  be  nearly  a  factor 
of  three  higher  where  the  interconnect  contacts  the  electrode 
than  midway  between  contacts,  regardless  of  the  axial 
position.  Resistance  to  radial  distribution  of  current  in  the 
LSM  cathode  was  found  to  be  very  high.  Effects  of  reduction 
of  tube  radius  or  changes  in  thickness  of  the  electrodes, 
however,  were  not  discussed.  In  general,  the  role  of  various 
parameters  of  the  interconnect  design  is  difficult  to  visualize 
with  numerical  models,  except  in  selected  cases. 

The  purpose  of  this  work  is  to  present  analytical  expres¬ 
sions,  which  describe  the  dependence  of  the  area- specific 
resistance  (ASR)  of  a  SOFC  stack-repeat  unit,  consisting 
of  a  cell  and  an  interconnect,  as  a  function  of  the  cathode 
and  anode  thicknesses,  charge  transfer  resistances  (activa¬ 
tion  polarization)  at  the  cathode-electrolyte  and  anode- 
electrolyte  interfaces,  ionic  resistivity  of  the  electrolyte, 
electronic  resistivities  of  the  anode  and  the  cathode,  respec¬ 
tively,  distance  between  interconnect  contacts,  and  inter¬ 
connect-electrode  contact  area.  The  analysis  treats  both 
one-  and  two-dimensional  interconnect  symmetries.  The 
area-specific  resistance  (ASR)  will  be  used  as  the  figure  of 
merit  to  assess  the  design  of  a  repeat  unit  consisting  of  a 
single-cell  and  interconnect.  It  will  not  be  the  objective  of 
this  paper  to  include  temperature  variations  or  the  effects  of 
fuel  utilization.  This  will  be  the  subject  of  a  future  inves¬ 
tigation. 


2.  Analysis 

There  are  two  basic  interconnect  symmetries,  which  are 
used  in  planar  SOFC  stacks,  one-dimensional  (channels)  and 
two-dimensional  (dimples).  Fig.  1  shows  an  SOFC  cell- 
interconnect  repeat  unit  with  a  one-dimensional  intercon¬ 
nect  symmetry.  The  interconnect  contacts  of  width  2v0  are 
periodically  spaced  at  a  distance  of  L,  and  their  length 
extends  across  the  entire  cell.  Similarly,  Fig.  2  shows  another 
SOFC  cell-interconnect  repeat  unit  with  a  two-dimensional 
interconnect  symmetry.  The  interconnect  contacts  are  cir¬ 
cular  in  shape  and  are  placed  in  a  hexagonal  array  to  cover 
the  surface  of  the  electrode.  The  center-to-center  distance 
between  nearest  neighbor  interconnect  contacts,  as  shown  in 
the  figure,  is  M;  and  the  contact  area  per  contact  with  the 
electrode  is  nr q. 

The  main  objective  of  this  analysis  is  to  calculate  the 
area-specific  resistance  per  repeat  unit  consisting  of  a  cell 
and  interconnect  for  both  interconnect  designs.  Once  this  is 
obtained,  the  stack  resistance  is  determined  by  simply 
multiplying  the  area-specific  resistance  of  the  repeat  unit 
by  the  number  of  cells  in  the  series-connected  stack  and 
dividing  by  the  active  area  per  cell.  It  is  assumed  at  the 
outset  that  all  transport  properties  are  ohmic.  Non-ohmic 
transport  can  be  treated,  but  this  would  require  numerical 
solution  of  the  steady-state  transport  equations.  It  is  also 
assumed  that  transport  of  gaseous  species  to  and  from  the 
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Fig.  2.  A  schematic  of  a  SOFC  in  series  contact  with  interconnect  of  a  two-dimensional  geometry. 


electrodes,  as  well  as  through  the  electrodes,  is  fast.  That  is, 
concentration  polarization  is  assumed  to  be  negligible. 
These  assumptions  are  justifiable  for  two  reasons:  (1) 
the  current  versus  voltage  response  of  SOFCs  (in  stacks) 
is  usually  nearly  linear,  and  (2)  it  has  been  shown,  both 
theoretically  and  experimentally,  that  concentration  polar¬ 
ization  becomes  significant  only  for  thick  electrodes 
(>1  mm)  operating  at  high  current  densities. 

In  a  steady-state,  the  divergence  of  the  flux  of  each  species 
is  zero,  and  for  spatially  invariant,  ohmic  transport  proper¬ 
ties,  this  leads  to  the  result  that  the  electrostatic  potential 
within  all  regions  of  a  SOFC  must  satisfy  the  Laplace 
equation,  V2</>  =  0.1  In  principle,  the  Laplace  equation 
should  be  solved  subject  to  appropriate  boundary  conditions 
within  each  region  of  the  SOFC  (cathode,  electrolyte,  and 
anode).  The  complexity  of  the  boundary  conditions  and 
number  of  regions  of  various  materials,  make  analytic 
expressions  for  the  area-specific  resistance  based  upon 
solutions  to  Laplace  equations  difficult.  Analytic  expres¬ 
sions,  however,  can  be  obtained  using  a  ladder  network 
approach  if  two  assumptions  are  made.  (1)  The  electrostatic 
potential  within  the  electrocatalyst  materials  of  the  cathode 


1  This  assumes  that  each  region  is  locally  homogeneous  and  a 

continuum,  and  that  electrochemical  reactions  are  confined  to  thin  regions 
at  electrolyte/electrode  interfaces. 


and  anode  is  invariant  in  the  vertical  direction,  namely,  dcj)/ 
dy  =  0.  This  is  justified  if  the  resistivities  of  the  electrodes 
are  much  smaller  than  that  of  the  electrolyte.  In  a  SOFC,  the 
materials  used  for  the  electrocatalysts  in  the  cathode  (LSM, 
pc  =  0.01  Q  cm  and  anode  pa  =  0.001  Q  cm)  have  resistiv¬ 
ities  that  are  significantly  lower  than  that  of  the  electrolyte 
(eight  YSZ,  pe  =  50  Q  cm).  (2)  In  the  electrolyte,  the  elec¬ 
trostatic  potential  exhibits  a  much  greater  variation  through 
the  thickness  rather  than  along  the  plane  of  the  electrolyte, 
i.e.  d(j)/dy>  d(j)/dx.  Since  the  resistivities  of  the  electro¬ 
catalyst  materials  in  a  SOFC  are  rather  small,  the  electro¬ 
static  potentials  in  the  cathode  and  the  anode  vary  slowly  in 
the  horizontal  (in-plane)  direction.  Since  the  electrolyte  lies 
between  the  cathode  and  the  anode,  the  electrostatic  poten¬ 
tial  within  the  electrolyte  varies  slowly  in  the  horizontal 
direction,  as  well.  This  implies  that  current  flows  from  the 
anode  to  the  cathode  through  the  electrolyte  primarily  in  the 
vertical  (through  the  thickness)  direction. 

2.7.  One -dimensional  interconnect  symmetry 

The  one-dimensional  stack-repeat  unit  used  for  the  ana¬ 
lysis  is  shown  in  Fig.  3.  The  potentials  at  the  interconnect- 
electrode  interfaces  are  fixed.  Currents  are  shown  entering 
and  leaving  differential  elements  within  the  cathode  and  the 
anode.  In  steady- state,  currents  entering  and  leaving  any 
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x=0 


Fig.  3.  A  schematic  showing  a  cross-section  of  the  SOFC  depicted  in  Fig.  1.  This  view  is  obtained  by  considering  the  symmetry  associated  with  the 
periodically  repeated  one-dimensional  interconnect  geometry. 


differential  element  must  be  conserved.  This  condition  for 
an  element  in  the  cathode  gives 


tc 

Pc 


4> i  {x  +  Ax)  -  (j) i  (x)~ 

Ax 

Ax 

(~E 

.Pe^  +  Rct  + 

k 

Pc  L 


</>i  (x  —  Ax)  —  </>j  (x) 


Ax 


(~E  ~  4>i(x)  +  <t>i(x))  =  0  (1) 


The  parameters  in  this  equation  are:  the  position  measured 
horizontally  from  the  center  of  the  interconnect  contact,  x; 
the  electrostatic  potential  in  the  cathode  as  a  function  of 
position,  </>j  (x) ;  the  electrostatic  potential  in  the  anode  as  a 
function  of  position,  </)j  (x);  the  thickness  of  the  cathode,  4; 
the  thickness  of  the  electrolyte,  4;  the  thickness  of  the  anode 
4;  the  electronic  resistivity  of  the  cathode  (assumed  for  the 
moment  to  be  constant  through  the  thickness),  pc;  the  ionic 
resistivity  of  the  electrolyte,  pe;  the  electronic  resistivity  of 
the  anode  (assumed  to  be  constant  through  the  thickness), 
pa;  the  charge  transfer  resistance  at  the  cathode-electrolyte 
interface,  Rcct ;  the  charge  transfer  resistance  at  the  anode- 
electrolyte  interface,  7^t;  and  the  Nernst  potential,  E. 
Although  for  the  purposes  of  simplifying  the  analysis,  the 


charge  transfer  resistances  are  given  as  being  associated  with 
the  physically  distinct  electrode-electrolyte  interfaces,  it 
can  be  shown  that  the  charge  transfer  resistances  themselves 
are  functions  of  the  electrode  structure  [3-6].  This  aspect 
will  be  elaborated  upon  later.  Taking  the  limit  as  Ax  —>  0  of 
the  first  two  terms  in  Eq.  (1)  and  then  dividing  by  Ax  gives 


4 

"d  4>l(x) 

d 

-E  -  $(x)  +  <fi(x) 

pcAx 

dx 

+ 

dx 

__ 

Pete  +  +  Rct 

(2) 

Again,  taking  the  limit  as  Ax  — >  0  and  rearranging  terms 
gives 

-^^-+ko[-E-m+^(x)]=o  o) 

where 


kc  = 


PcAc 


Pje  +  R%  +  R% 


(4) 


If  the  cathode  is  graded  such  that  the  electronic  conductivity 
is  a  function  of  through  the  thickness,  the  pc  from  the  above 
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equations  is  replaced  by  the  cathode  resistivity,  (pc),  given  boundary  conditions  applicable  to  the  geometry  of  Fig.  3 
by  are 


(Pc) 


tC 

fody/pdy) 


(5) 


where  pc(y )  is  the  electronic  resistivity  of  the  cathode  at  a 
position  y,  measured  with  respect  to  one  of  the  interfaces; 
cathode/electrolyte  interface,  or  cathode/oxidant  gas  phase 
interface.  For  a  graded  anode,  similar  average  value  of 
anode  resistivity,  (pa),  can  be  defined.  Differential  Eq.  (3) 
describes  the  spatial  dependence  of  the  electrostatic 
potential  within  the  cathode  in  the  range  xo  <  x  <  L/2  for 
the  one-dimensional  interconnect  symmetry.  There  are  two 
unknowns  in  Eq.  (3),  namely,  the  electrostatic  potentials 
within  the  cathode  and  the  anode,  as  functions  of  position.  In 
order  to  solve  Eq.  (3),  the  spatial  dependence  of  the  electro¬ 
static  potential  within  the  anode  must  also  be  evaluated. 
In  steady- state,  the  current  within  a  differential  element  in 
the  anode  must  also  be  conserved  (see  Fig.  3).  Repeating  the 
same  procedure,  as  was  used  for  the  cathode,  gives 

+  ka[E  -  #(x)  +  0,a(x)]  =  0  (6) 


where 


PaAa 

Pe4  + 


(7) 


In  the  case  of  anode  of  graded  properties,  pa  from  the  above 
equations  is  replaced  by  (pa)  which  is  defined  by  an  equation 
analogous  to  Eq.  (4).  Differential  Eqs.  (3)  and  (6)  are  linearly 
independent  and  can,  in  principle,  be  solved  simultaneously 
to  obtain  </>j  (x)  and  </>j  (x) .  In  order  to  reduce  this  system  of 
two  differential  equations  and  two  unknowns  to  a  single 
differential  equation  involving  only  one  unknown,  the  fol¬ 
lowing  substitution  will  be  made: 


viM  =  $(*)  -  (8) 

Eq.  (8)  describes  the  difference  in  electrostatic  potentials 
between  the  cathode  and  the  anode  at  a  given  position,  v. 
Subtracting  Eq.  (6)  from  Eq.  (3),  and  using  the  above 
substitution  gives 


1 

kc  +  K 


d2vi(x) 

dx2 


+  Vi(x)  =  E 


(9) 


Vi(x0)  =  </>° 


and 


dvi(x) 

dx 


=  0 

x=L/2 


(12) 


The  first  boundary  condition  corresponds  to  the  fixed  poten¬ 
tials  at  the  interconnect  contacts.  The  mirror  plane  at  v  =  L/2 
across  which  no  current  flows  gives  the  second  boundary 
condition  (symmetry).  Solving  for  V!  and  v2  gives 


vi  =  v2  e  aL  and  v2 


0°  -E 

ea(x0-L)  _|_  Q-ax0 


(13) 


Substitution  of  Eq.  (13)  into  Eq.  (10)  gives  the  complete 
solution  for  the  difference  in  electrostatic  potentials  between 
the  cathode  and  the  anode,  v:(x),  as  a  function  of  position 
when,  xo  <  x  <  L/2 ,  namely, 


Vi(x)  =  E  +  0° 


E) 


g a(x—L ) 

qo(xq—L)  q- ax o 


(14) 


This  equation  may  be  used  as  an  intermediate  result  to  obtain 
solutions  for  </>j(x)  and  (j>l(x).  Alternatively,  it  is  sufficient 
for  the  evaluation  of  the  area-specific  resistance. 

Examination  of  Fig.  3  shows  that  there  are  two  separate 
regions  of  the  electrolyte  that  contribute  to  the  current  flow, 
/j.  For  the  inner  region  located  directly  beneath  the  inter¬ 
connect  contact,  that  is  for  0  <  x  <  xq,  the  potential  differ¬ 
ence  is  constant  and  is  given  by  Eq.  (12).  From  Fig.  3,  the 
current  flowing  through  the  electrolyte  in  the  inner  region, 
per  unit  length  in  the  z-direction  (into  the  paper)  is 


(E  -  0°)xo 
/Ve  +  RCct  +  Ret 


(15) 


and  is  over  the  region  0  <  x  <  xq.  In  the  outer  region,  that  is, 
for  xq  <  x  <  L/2 ,  the  potential  difference  between  the  cath¬ 
ode  and  the  anode  varies  as  a  function  of  position.  From 
Fig.  3,  the  current  flowing  through  a  differential  element  of 
the  electrolyte  is  given  by 


dJo  (£  -  vi(*)) 

1  Pe?e  +  RCa  +  R*ct 


(16) 


per  unit  length  in  the  z-direction  (into  the  paper).  This 
current  is  found  by  substituting  Eq.  (14)  into  Eq.  (16)  and 
integrating  from  v0  to  L/2,  and  is  given  by 


The  solution  to  this  differential  equation  is  easily  shown  to 
be 

vi(x)  =  E  +  vi  eax  +  v2  e~ax  (10) 

where 


/,= 


E  —  (j>° 


a(pete  +  RLa  +  R'a) 


Q-ax o  _  ga(xo-L)' 


Q-ax0  _|_  ea(x0-L) 


(17) 


The  net  current  (over  the  distance  0  to  L)  is  given  by  the  sum 
of  Eqs.  (15)  and  (17),  that  is, 


h=A  +  7I0  (18) 

Evaluation  of  the  area-specific  resistance  is  now  straightfor¬ 
ward.  The  resistance  of  the  structure  shown  in  Fig.  3 
and  Vi  and  v2  are  constants,  which  must  be  determined  between  the  interconnect  contacts  can  be  found  by  applying 

by  applying  the  requisite  boundary  conditions.  The  two  Ohm’s  law  to  the  expression  for  the  total  current,  i.e. 


a  = 


/  PcAc+PaAa  _  ,/j— nr 

V  Pete  +  Rcct  +  Ret  Vc  3 


(ID 
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absolute  value  of  ( d(j)°/dli ).  Interconnect  material  acts  in 
series  and  contributes  approximately  a  resistance  of  pfi/x o 
to  the  total.2  Multiplying  the  total  resistance  by  the  area  (per 
unit  dimension  along  the  channels),  LI 2,  gives  an  analytical 
expression  for  the  area-specific  resistance  (ASR)  of  the 
stack-repeat  unit,  given  by 


where  kc  is  the  same  parameter  as  defined  earlier.  Taking  the 
limit  as  Ar  — >  0  gives  the  differential  equation 


d20n(r) 


dr2 


dr 


+  rkc(—E  -  0n(r)  +  0„(r))  =  0 

(23) 


R  _L  f _ pe4  +  RLcl  +  R'^ _ Pjtj\ 

1  2  \x0  +  (l/a)[(e_ax°— ea(-lo_z'))/(e_a,:o+ea(xo_i))]i~  x0  / 

(19) 

in  units  of  Q  cm2.  Note  that  this  equation  is  a  function  of  the 
various  material  and  design  parameters. 

2.2.  Two-dimensional  interconnect  symmetry 


The  analysis  of  the  two-dimensional  interconnect  sym¬ 
metry  will  proceed  in  a  manner  similar  to  the  one  presented 
here  for  the  one-dimensional  symmetry.  The  repeat  unit  (in 
the  plane  of  the  interconnect)  for  the  two-dimensional 
interconnect  symmetry,  shown  in  Fig.  2,  is  a  hexagon  with 
sides  of  length  M/2.  For  the  purpose  of  simplifying  the 
mathematics,  the  hexagonal  repeat  unit  will  be  approxi¬ 
mated  by  a  circle  which  introduces  a  cylindrical  symmetry. 
The  radius  of  this  circle,  Rl 2,  will  be  chosen  such  that  it  has 
an  area  equal  to  that  of  the  hexagon.  The  relationship 
between  R  and  M  is 


R  =M 


1/4 


(20) 


Fig.  4  shows  a  circular  section  of  a  SOFC  with  the  inter¬ 
connect  contact  located  at  the  center.  The  figure  shows 
current  entering  and  leaving  a  cylindrical  ring  of  thickness 
Ar  in  both  the  cathode  and  the  anode.  In  steady-state,  the 
current  entering  and  leaving  the  ring  within  the  cathode  must 
be  conserved.  This  leads  to 


Differential  Eq.  (23)  describes  the  spatial  dependence  of 
electrostatic  potential  within  the  cathode  for  the  two-dimen¬ 
sional  interconnect  symmetry  in  the  range  ro  <  r  <  Rl 2. 
This  equation  involves  the  electrostatic  potential  within  the 
anode  as  a  function  of  position,  which  is  not  known.  The 
governing  equation  for  the  anode  is  found  in  a  similar 
manner,  and  is  given  by 


d2</>n(r)  d$jj(r) 


dr2 


dr 


+  rhJ-E  -  0n(r)  +  0„(r))  =  0 

(24) 


Differential  Eqs.  (23)  and  (24)  are  linearly  independent  and 
can  be  solved  simultaneously  to  obtain  </>n(r)  and  4>^(a).  In 
order  to  reduce  this  system  of  two  differential  equations  and 
two  unknowns  to  a  single  differential  equation  involving 
only  one  unknown,  the  following  substitution  will  be  made 

vnM  =  <£n  (r)  -  <£ii(r)  (25) 

Eq.  (25)  is  the  difference  in  electrostatic  potentials  between 
the  cathode  and  the  anode  at  a  position,  r.  Subtracting 
Eq.  (24)  from  Eq.  (23),  and  making  the  above  substitution 
gives 

_  d  vn6)  _  ldv”(r)  +  (£c  +  kc)vu(r)  =  (kc  +  ka)E  (26) 

This  equation  differs  from  Eq.  (9)  by  virtue  of  the  presence 
of  the  term  (l/r)(dvn(r)/dr).  Eq.  (26)  can  be  identified  as  the 
modified  Bessel  equation  [10,11].  The  solution  to  this 
differential  equation  is  given  by  [10,11] 


2  ntc(r  +  Ar/2) 
Pc 

2  7itc(r  —  Ar/2) 


0n(r  +  Ar)  </>n(r) 


Ar 

0n(r  —  Ar)  —  0n(r) 


7i(r+Ar/2)  — 7i(r  —  Ar/2 ) 
Pe'e+tfg  t+^ct 


Ar 

2 


(  E  0n(r)  +  0n(r))  —  ^ 

(21) 


Expanding  Eq.  (21),  rearranging  terms,  and  factoring  out 
terms  in  Ar  gives 


r 

Ar 


0n(r  +  Ar)  </>n(r)  0  n(r)  ^n(r  ^r) 


Ar 


0n 6  +  Ar)  -  (j)Ln(r  -  Ar) 


2  Ar 


Ar 


+  rkc(-E-  0„(r) 


+  0n(r))  —  o 


(22) 


2  This  assumes  tt  much  larger  than  the  thickness  of  interconnect  away 
from  the  ribs. 


vn(r)  =  E  +  v3I0(ar)  +  v4K0(ar )  (27) 

where  I0(ar )  and  K0(ar )  are  the  modified  Bessel  functions  of 
zeroth-order  of  the  first  and  second  kind,  respectively,  and  v3 
and  v4  are  constants.  The  constants  v3  and  v4  can  be  found  by 
applying  the  appropriate  boundary  conditions,  namely, 


vii(t*o)  — 


and 


dvn(r) 

dr 


=  0 

r=R/ 2 


(28) 


Again,  the  first  boundary  condition  represents  the  fixed 
potential  between  the  interconnect  contacts.  The  second 
boundary  condition  is  obtained  by  deducing  that  no  current 
flows  across  the  perimeter  at  r  =  RI2.  It  can  be  shown  that 
the  constants  v3  and  v4  are  given  by 


v3 


4>°  -E 


v4  =  v3 


[Io(ar0)  +  {7i(aR/2)/^i(flK/2)}«b(an,)] 
I\  (aR/2) 


and 


K\  (aR/2) 


(29) 
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Fig.  4.  A  schematic  showing  a  section  of  the  SOFC  depicted  in  Fig.  2.  This  circular  section  is  obtained  by  approximating  the  hexagonal,  in-plane  repeat  unit 
of  the  two-dimensional  interconnect  geometry  with  a  circle  of  the  same  area. 


and  the  complete  solution  for  the  potential  drop  across  the 
cell  as  a  function  of  position  is  given  by 

v  fr)  F  |  /  .o  p\  (  +  Kfar)h(aR/2)  \ 

IlU  ’  \I0(aro)K\ [aR/ 2)  +  K0{ar0)h (aR/2) ) 

(30) 


For  the  inner  region  located  directly  beneath  the  intercon¬ 
nect  contact,  that  is  0  <  r  <  ro,  the  potential  difference  is 
constant,  and  the  current  flowing  is  given  by 


(£  -  <jP)nrl 

Pete  +  RLcl  +  /?c  t 


(31) 


In  the  above  equations,  K\  and  I\  are  respectively  modified 
Bessel  functions  of  first-order,  and  of  first  and  second  kind. 
Eq.  (30)  may  be  used  as  an  intermediate  result  to  obtain 
solutions  for  (j)\Y  (r)  and  (r).  Alternatively,  it  is  sufficient  for 
the  evaluation  of  the  area-specific  resistance  of  the  repeat  unit. 

Examination  of  Fig.  4  shows  that  there  are  two  separate 
regions  of  the  cell  that  contribute  to  the  total  current  flow,  fu. 


In  the  outer  region,  that  is  ro  <  r  <  Rt 2  the  difference  in 
potentials  between  the  cathode  and  the  anode  varies  as  a 
function  of  position.  The  current  flow  through  a  ring  of  the 
electrolyte  is  given  by 


<{r) 


2nr(E  —  vn(r))  dr 

/Ve  +  Rcc  t  +  Ret 


(32) 
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Substitution  for  vn(r)  from  Eq.  (30)  in  the  above  equation, 
and  integration  between  r0  and  R/2  gives  the  corresponding 
current,  given  by 

jO  _  (  27iro(£  -  £)  \ 

”  W/Vc  +  +  Ra)J 

Ki  (ar^Ii  (aR/2)  -  h  {ar^Ki  (gR/ 2)\ 
K0(ar0)h(aR/2)  +  I0(ar*)Ki(aR/2)J  1  ’ 

The  net  current  flowing  through  the  electrolyte  across  area  of 
tiR2IA  is  given  by 

In  =  I[ i  +  /',!  (34) 

The  resistance  of  the  structure  shown  in  Fig.  4,  for  a  pair  of 
opposite  interconnect  contacts,  and  the  corresponding  area, 
is  given  by  applying  Ohm’s  law,  i.e.  d</>°/d /n.  The  inter¬ 
connect  material  contributes  a  resistance  of  pfi / (nr2}),  to  the 
total.3  Multiplying  the  total  resistance  by  the  area,  7iR2/4, 
gives  an  analytical  expression  for  the  area-specific  resistance 
(ASR)  of  the  stack  per  one  repeat  unit,  and  is  given  by 


pm,  and  R*  =s-  0.1  Q  cm2.  All  transport  parameters  are  cho¬ 
sen  to  be  the  same  for  the  two  types  of  cells.  The  charge 
transfer  resistances  depend  upon  the  intrinsic  charge 
transfer  characteristics  (such  as  that  occurring  at  a  three 
phase  boundary,  TPB)  and  the  micro  structure  (e.g.  the 
TPB  length).  The  calculations  given  here  assume  that  the 
conductivity  of  the  interconnect  material  is  much  greater 
than  any  component  of  the  cell.  If  this  is  not  the  case,  the 
term  corresponding  to  the  resistance  of  the  interconnect 
material  must  be  included. 

Figs.  5  and  6  are  plots  of  the  calculated  area-specific 
resistance  for  one-  and  two-dimensional  interconnect  sym¬ 
metries  as  a  function  of  the  spacing  between  contacts,  L  or  R , 
for  the  electrolyte-  and  anode- supported  cells,  respectively. 
It  can  be  seen  from  the  plots  that  regardless  of  the  inter¬ 
connect  contact  spacing,  the  one-dimensional  symmetry 
always  gives  a  lower  area- specific  resistance  than  the 
two-dimensional  symmetry,  consistent  with  expectations. 
For  closely  spaced  interconnect  contacts,  the  area- specific 
resistance  is  approximately  equal  to  the  single-cell  resis- 


Rn 


^  f _ /Ve  +  RCct  +  Kt _ .  Pjtj\ 

4  Vo  +  (2r0/a)[{K1(ar0)I1(aR/2)  -  h{arQ)Kx{aR/2)} / {K^arMaR/2)  +  /0(flr0)tf1(aR/2)}]  r*  ) 


(35) 


Note  that  this  is  a  function  of  the  various  material  and  design 
parameters. 

3.  Results  and  discussion 

Eqs.  (19)  and  (35)  give  the  area- specific  resistances  (ASR) 
in  terms  of  transport  and  design  parameters  of  a  SOFC  stack- 
repeat  unit  for  one-  and  two-dimensional  interconnect  sym¬ 
metries.  They  show  that  it  is  not  possible  to  design  an 
interconnect  independent  of  cell  parameters  since  the  overall 
area-specific  resistance  depends  upon  a  number  of  para¬ 
meters  pertaining  to  the  cell.  Two  types  of  planar  cell 
designs  are  currently  under  investigation  at  several  labora¬ 
tories:  (1)  thick  YSZ  electrolyte  as  the  support  member  with 
thin  electrodes,  screen-printed  onto  the  electrolyte,  and  (2) 
thin  YSZ  electrolyte- supported  on  a  thick  porous  anode  with 
a  cathode  of  intermediate  thickness.  In  what  follows,  calcu¬ 
lations  are  presented  for  typical  values  of  various  parameters 
for  both  types  of  cells.  We  will  refer  to  these  as  an  electrolyte- 
supported  cell  and  an  anode-supported  cell,  respectively.  The 
parameters  chosen  for  the  electrolyte-supported  cell  are: 
pe  =  50  Q  cm,  4  =  150  pm,  pc  =  0.1  Q  cm,  tc  =  25  pm, 
Rcct  =  0.2  Q  cm2,  pa  =  0.001  Q  cm,  4  =  25  pm,  and  R^t  = 
0.1  Q  cm2.  The  interconnect  contact  dimension  is  selected  as 
v0  or  r0  =  10  pm.  The  parameters  chosen  for  the  anode- 
supported  cell  are:  pe  =  50  Q  cm,  4  =  10  pm,  pc  = 
0.1  Q  cm,  4  =  100  pm,  Rcct  =  0.2  Q  cm2,  pa  =  0.001  Q  cm, 


3  This  assumes  that  tt  is  much  greater  than  interconnect  thickness  away 
from  the  dimples. 


tance,  pe4  +  Rcct  +  R^v  As  the  interconnect  contact  spacing 
is  increased,  a  transition  in  behavior  occurs.  For  large 
interconnect  contact  spacing,  the  area-specific  resistance 
increases  due  to  the  resistance  to  flow  of  electrons  in  the 


Electrolyte-Supported  Cell 


Fig.  5.  Plots  of  the  area-specific  resistance  (electrolyte-supported)  for  one- 
and  two-dimensional  interconnect  geometries,  respectively,  as  a  function 
of  spacing  between  the  interconnect  contacts,  L  or  R.  The  values  used  in 
the  calculation  for  an  electrolyte- supported  cell  are:  pe  =  50  Q  cm,  te  = 
15  pm,  pc  =  0.1  Q  cm,  tc  =  25  pm,  Rcct  =  0.2  f2  cm2,  pa  =  0.001  Q  cm, 
ta  =  25  pm,  R^t  =  0.1  0  cm2,  and  x0  or  r0  =  10  pm. 
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Anode -Supported  Cell 


Fig.  6.  Plots  of  the  area-specific  resistance  (anode-supported)  for  one-  and 
two-dimensional  interconnect  geometries,  respectively,  as  a  function  of 
spacing  between  the  interconnect  contacts,  L  or  R.  The  values  used  in  the 
calculation  for  an  anode- supported  cell  are:  pe  =  50  Cl  cm,  te  —  10  pm, 
pc=0.1Qcm,  tc  =  100  pm,  =  0.2  Q  cm2,  pa  =  0.001  Q  cm, 

4  =  500  pm,  R^t  =  0.1  Q  cm2,  and  x0  or  r0  =  10  pm. 

electrocatalyst  materials  within  the  electrodes,  i.e.  electrode 
sheet  resistance.  The  transition  from  single-cell-dominated 
to  interconnect  design-dominated  resistance  occurs  at  a 
lower  interconnect  contact  spacing  for  the  two-dimensional 
symmetry.  The  plot  of  Eq.  (19)  for  the  one-dimensional 
symmetry  shows  that  for  widely  spaced  interconnect  con¬ 
tacts,  the  area-specific  resistance  increases  linearly  with 
interconnect  contact  spacing,  L.  By  comparison,  the  plot 
of  Eq.  (35)  for  the  two-dimensional  symmetry  shows  that  for 
widely  spaced  interconnect  contacts,  the  area- specific  resis¬ 
tance  increases  parabolically  with  interconnect  contact  spa¬ 
cing,  R.  Also,  interconnect  contact  spacing  for  designs  with 
a  two-dimensional  symmetry  must  be  smaller  than  for 
designs  with  a  one-dimensional  spacing  to  achieve  the  same 
area-specific  resistance.  As  a  numerical  example,  if  the 
desired  maximum  repeat  unit  ASR  is  0.5  Q  cm2  for 
anode- supported  cell  design  (Fig.  6),  the  interconnect  con¬ 
tact  spacing  for  two-dimensional  geometry  must  be  under 
2  mm  (<0.2  cm),  while  that  for  the  one-dimensional  geo¬ 
metry  can  be  as  large  as  ~5  mm  (~0.5  cm). 

Comparison  of  Figs.  5  and  6  also  shows  that  for  the  same 
transport  properties,  the  anode- supported  cell  exhibits  a 
lower  area- specific  resistance  than  the  electrolyte- supported 
cell.  This,  of  course,  is  to  be  expected,  as  the  electrolyte 
resistance  is  higher  in  the  case  of  electrolyte- supported  cells. 

In  the  calculations  given  in  Figs.  5  and  6,  the  cathode 
resistivity  was  assumed  to  be  0.1  Q  cm.  The  typical  value  of 
FSM  is  about  0.01  Q  cm.  However,  as  the  cathode  is  porous, 


Electrolyte-Supported  Cell 


Fig.  7.  Plots  of  the  area-specific  resistance  (electrolyte-supported)  for  one- 
and  two-dimensional  interconnect  geometries,  respectively,  as  a  function 
of  cathode  resistivity,  pc.  The  values  used  in  the  calculation  for  an 
electrolyte-supported  cell  are:  pe  =  50  Q  cm,  te  =  150  pm,  tc  =  25  pm, 
Rcct  =  0.2  Q  cm2,  pa  =  0.001  Q  cm,  ta  =  25  pm,  R^t  =  0.1  Q  cm2,  x0  or 
ro  =  10  pm,  and  L  or  R  =  3  mm. 

generally  contains  some  YSZ,  and  that  its  conductivity  can  be 
lower  in  oxygen-depleted  air,  such  as  during  the  actual  SOFC 
operation,  the  actual  cathode  resistivity  may  be  much  higher 
than  0.01  Q  cm,  hence,  the  choice  of  0.1  Q  cm.  Further,  if 
FSM  reacts  with  YSZ  to  form  Fa2Zr207,  the  resistivity  may 
be  even  higher  than  0.1  Q  cm.  However,  if  the  cathode 
consists  of  materials  such  as  cobaltites,  and  their  reaction 
with  YSZ  can  be  suppressed  (by  using  a  ceria  buffer  layer)  the 
resistivity  could  be  much  lower  than  0.01  Q  cm.  Thus,  for  the 
purposes  of  evaluating  the  role  of  cathode  resistivity,  calcula¬ 
tions  of  the  area-specific  resistance  (ASR)  were  made  over  a 
range  of  cathode  resistivities,  for  both  electrolyte-supported 
and  anode-supported  cells,  and  for  both  interconnect  geome¬ 
tries,  for  an  interconnect  contact  spacing  of  3  mm,  and 
interconnect  contact  dimension  of  10  pm.  The  results  of  these 
calculations  are  shown  in  Figs.  7  and  8.  Note  that  as  the 
cathode  resistivity  rises,  the  ASR  rises;  and  once  again,  it  rises 
faster  for  the  two-dimensional  symmetry.  For  the  two-dimen¬ 
sional  geometry,  in  particular,  the  rise  in  ASR  is  rather  rapid, 
emphasizing  the  importance  of  the  need  for  a  low  sheet 
resistance.  For  the  anode-supported  design  (Fig.  8),  which 
has  lower  overall  ASR  compared  to  the  electrolyte-supported 
design,  note  that  the  repeat  unit  ASR  for  the  two-dimensional 
geometry  increases  from  ~0.4  Q  cm2  for  pc  =  0.01  Q  cm,  to 
~0.8  Q  cm2  for  pc  =  0.1  Q  cm.  That  is,  the  stack  perfor¬ 
mance  drops  by  a  factor  of  2.  Beyond  pc  =  0. 1  Q  cm,  the  rise 
is  even  more  rapid.  Over  the  same  range,  the  one-dimensional 
geometry  is  more  forgiving. 
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Anode-Supported  Cell 


Fig.  8.  Plots  of  the  area-specific  resistance  (anode-supported)  for  one-  and 
two-dimensional  interconnect  geometries,  respectively,  as  a  function  of 
cathode  resistivity,  pc.  The  values  used  in  the  calculation  for  an  anode- 
supported  cell  are:  pe  =  50  £2  cm,  4  =  10  pm,  4  =  100  pm, 
Rcct  =  0.2  Q  cm2,  pa  =  0.001  £2  cm,  4  —  500  pm,  R*t  =  0.1  Q  cm2,  x0  or 
ro  =  10  pm,  and  L  or  R  =  3  mm. 


Electrolyte-Supported  Cell 


Interconnect  Contact  Half  Width  or  Radius  (pm) 


Fig.  9.  Plots  of  the  area-specific  resistance  (electrolyte- supported  cell)  for 
one-  and  two-dimensional  interconnect  geometries,  respectively,  as  a 
function  of  interconnect  contact  half-width  or  radius,  x0  or  r0.  The 
values  used  in  the  calculation  for  an  electrolyte- supported  cell  are: 
pe  =  50  Q  cm,  4  =  150  pm,  pc  =0.1  Qcm,  4  =  25  pm,  =  0.2  Q 
cm2,  pa  =  0.001  £2  cm,  4  =  25  pm,  R^t  =  0.1  Q  cm2,  L  or  R  =  3  mm. 


Figs.  9  and  10  are  plots  of  the  calculated  area-specific 
resistance  for  one-  and  two-dimensional  symmetries  as  a 
function  of  the  interconnect  contact  half- width  or  radius,  x0 
or  r0,  for  stacks  with  electrolyte- supported  and  anode- 
supported  cells,  respectively,  for  interconnect  contact  spa¬ 
cing  of  3  mm.  Once  again,  it  can  be  seen  from  the  figures 
that  the  area-specific  resistance  for  the  interconnect  contact 
design  with  a  one-dimensional  symmetry  is  always  lower 
than  the  design  with  a  two-dimensional  symmetry.  For 
large  interconnect  contact  areas,  the  area-specific  resis¬ 
tance  is  approximately  equal  to  the  single-cell  resistance,4 
pefe  +  Rcct  +  R4t,  and  as  the  contact  area  decreases,  the  area- 
specific  resistance  increases.  The  area- specific  resistance 
for  one-dimensional  symmetry  increases  to  an  asymptotic 
maximum  as  the  contact  half- width  (area)  decreases.  An 
examination  of  Eq.  (19)  shows  that  for  a  highly  conducting 
interconnect  material  (p,  —>  0),  and  small  rib  thickness 
(h  — >  0),  the  second  term  in  parenthesis  can  be  neglected. 
In  all  calculations  presented  here,  this  assumption  was 
made.  If  now  one  examines  the  limit  xq  — >  0,  the  repeat 
unit  ASR  becomes 


Anode-Supported  Cell 


aL(pete  +  Rcct  +  Rlt) 
2[(1  -  e-flL)/(l  +e-flL)] 


(36) 


4  Of  course,  interconnect  contact  area  cannot  be  increased  indefinitely, 
otherwise  there  would  be  no  space  left  for  the  transport  of  gaseous  species 
along  the  interconnect  faces. 


Interconnect  Contact  Half  Width  or  Radius  (pm) 

Fig.  10.  Plots  of  the  area-specific  resistance  (anode- supported  cell)  for 
one-  and  two-dimensional  interconnect  geometries,  respectively,  as  a 
function  of  interconnect  contact  half-width  or  radius,  x0  or  r0.  The 
values  used  in  the  calculation  for  an  anode- supported  cell  are:  pe  =  50  £1 
cm,  4  =  10  pm,  pc  =0.10  cm,  4  =  100  pm,  Rcct  =  0.2  £2  cm2,  pa  = 
0.001  £1  cm,  4  =  500  pm,  R^t  =  0.1  £1  cm2,  L  or  R  =  3  mm. 
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That  is,  the  ASR  approaches  a  finite  value.  A  similar 
analytical  expression  of  the  limit  of  Eq.  (35)  as  ro  — >  0,  is 
difficult  to  obtain.  It  can  be  shown,  however,  that  as  r0  — >  0, 
both  K0(ar0)  and  Ki(ar0 )  approach  infinity,  with  Ki(ar0) 
rising  faster.  If  the  ratio  Ki(aro)/K0(ar0 )  rises  slower  than  1/ 
r0,  one  would  expect  the  limit  of  Eq.  (35)  to  approach  infinity 
as  ro  — >  0.  Numerical  calculations  given  in  Figs.  9  and  10 
show  at  the  ASR  monotonically  increases  as  r0  decreases. 
That  is,  the  area-specific  resistance  for  a  two-dimensional 
symmetry  does  not  increase  to  a  finite  maximum  as  the 
contact  area  decreases.  In  fact,  Eq.  (35)  for  the  two-dimen¬ 
sional  area-specific  resistance  is  singular,  i.e.  infinite,  in  the 
limit  as  the  contact  radius,  r0  approaches  zero.  Thus,  the 
interconnect  contact  design  with  a  two-dimensional  sym¬ 
metry  can  have  a  significantly  higher  area- specific  resistance 
than  the  design  with  a  one-dimensional  symmetry  in  cases 
where  the  contact  area  is  small  and  the  spacing  between 
contacts  is  large.  Generally,  however,  the  interconnect  con¬ 
tact  dimension  is  expected  to  be  larger  than  10  pm.  Thus,  the 
interconnect  contact  dimension  is  not  expected  to  exhibit  a 
large  effect.  However,  in  an  actual  SOFC  stack,  due  to 
dimensional  irregularities,  it  can  often  so  happen  that  the 
contact  between  cell  and  interconnect  may  not  be  perfect. 
This  may  lead  to  a  situation  where  in  some  of  the  regions, 
interconnect  may  not  be  in  contact  with  the  cell.  This  would 
be  equivalent  to  locally  increasing  L  or  R ,  and  thus,  effec¬ 
tively  increasing  the  ASR.  The  analysis  presented  here 
shows  that  the  adverse  effect  of  irregular  contact  in  general 
would  be  greater  in  the  case  of  two-dimensional  geometry 
than  in  the  case  of  one-dimensional  geometry. 

The  interconnect  contact  spacing  at  which  area-specific 
resistance  (ASR)  transitions  from  single-cell  dominated  to 
sheet  (electrode)  resistance  dominated  can  be  estimated 
from  Eq.  (11).  For  interconnect  spacing  wider  than  ~l/a, 
resistance  to  distribution  of  current  in  the  electrodes  begins 
to  become  larger  than  the  resistance  of  the  single-cell, 
pe4  +  Rcc t  +  R^v  In  fact,  the  parameter  1  /a  arises  indepen¬ 
dent  of  the  boundary  conditions  or  shape  of  the  single-cells. 
As  a  general  rule  or  a  design  guideline,  the  spacing  between 
the  interconnect  contacts  should  be  less  than  1  /a  so  that 
performance  of  the  stack  reflects  performance  of  the  single¬ 
cells.  The  Siemens-Westinghouse  fuel  cell  design  modeled 
by  Bessette  et  al.  can  be  approximately  described  in  terms  of 
the  one-dimensional  interconnect  geometry  and  analyzed  in 
terms  of  the  parameter  1  /a  [9].  An  analytical  expression  for 
the  area- specific  resistance  similar  to  Eq.  (19)  can  be 
obtained  by  applying  slightly  different  boundary  conditions 
to  the  general  solutions  for  potential  as  a  function  of  posi¬ 
tion.  The  cell  parameters  used  by  Bessette  et  al.  are  as 
follows:  pe  =  9.98  Q  cm,  tG  =  40  pm,  pc  =  0.013  Q  cm, 
tc  =  7000  pm  (0.7  cm),  Rcct  +  R^t  =  0.7  Q  cm2  (estimate), 
pa  =  0.001  Q  cm,  and  4  =  100  pm.  For  these  parameters, 
the  value  of  1  /a  is  2.42  cm.  The  diameter  of  the  tube  on 
which  the  single-cell  is  built  was  assumed  to  be  1 1  cm  in  the 
work  of  Bessette  et  al.  [9].  Of  the  circumference,  ^28.4  cm 
is  active  electrolyte,  and  the  repeat  unit  length  consistent 


with  the  one-dimensional  interconnect  symmetry  is  half 
this  value,  14.2  cm.  Thus,  the  internal  resistance  of  the 
tube-based  design  assumed  in  the  calculations  by  Bessette 
et  al.  is  dominated  by  sheet  resistance  in  the  electrodes  [9]. 
Reduction  of  the  diameter  of  the  tubes  to  less  than  1.5  cm 
and  increase  in  thickness  of  the  electrode  should  give 
significantly  lower  area- specific  resistance. 

For  the  calculations  presented  in  Figs.  5-8,  the  chosen 
charge  transfer  resistances  correspond  to  composite  electro¬ 
des  that  are  mixtures  of  an  electrocatalyst  and  an  oxygen  ion 
conductor.  For  the  cathode,  for  example,  a  suitable  mixture 
is  that  of  Sr-doped  FaMn03  (FSM)  and  YSZ.  For  the  anode, 
a  suitable  mixture  is  Ni  +  YSZ.  A  porous  composite 
electrode  that  consists  of  a  contiguous  mixture  of  an  elec¬ 
trocatalyst  and  an  ionic  conductor  exhibits  a  much  lower 
overall  charge  transfer  resistance  than  a  single-phase,  elec¬ 
tronically  conducting  electrocatalyst  (e.g.  FSM).  This 
improved  performance  of  a  composite  electrolyte  is  largely 
due  to  increased  TPB  length,  which  leads  to  a  spreading  of 
the  reaction  zone  into  the  electrode.  Distance  to  which  this 
spreading  occurs  depends  upon  a  number  of  parameters, 
including  the  micro  structure.  In  general,  the  finer  the  micro¬ 
structure,  the  lower  is  the  overall  charge  transfer  resistance, 
the  better  is  the  performance,  and  the  smaller  is  the  required 
thickness  of  the  electrode  (from  the  standpoint  of  lower 
polarization). 

Recent  work  has  shown  that  in  terms  of  the  intrinsic 
charge  transfer  resistance  of  the  cathode,  Rcco0y  which  for 
example  describes  electrocatalytic  behavior  of  FSM  (of  a 
given  FSM  particle  size  and  particle  number  density  per  unit 
area  of  YSZ)  on  YSZ,  the  effective  charge  transfer  resistance 
of  a  composite  cathode  comprising  a  porous  mixture  of  FSM 
and  YSZ,  namely,  RccV  is  given  by  [3,6] 


dcPlK 


et(0) 


1  ~  Pc 


(37) 


where  pc  is  the  cathode  porosity,  pce  is  the  ionic  resistivity  of 
the  composite  cathode  (that  of  the  YSZ  in  the  composite 
cathode),  and  dc  is  the  particle  size  of  the  YSZ  in  the 
composite  cathode.  The  above  equation  is  applicable  typi¬ 
cally  for  cathode  thicknesses  in  excess  of  a  critical  thickness, 
2,  defined  by  [6] 


2  = 


^c^ct(O)  ( 1  -Pc) 


(38) 


More  specifically,  Eq.  (37)  is  applicable  if  the  electrode 
thickness  is  greater  than  about  two  or  three  times  2.  The 
typical  grain  size  of  YSZ  in  the  FSM-YSZ  cathode  is  dc  ~  1 
micron.  For  a  typical  volume  fraction  of  FSM  (~0.5)  and 
particle  size  (~1  pm),  an  approximate  value  of  ^t(o)  at 
800  °C  is  about  1  Q  cm2.  Finally,  for  YSZ  in  the  cathode, 
the  ionic  resistivity,  pce  is  about  50  Q  cm.  For  a  cathode 
porosity  of  ~0.3,  using  the  values  mentioned  here,  Eq.  (38) 
gives  2  «  12  microns.  Thus,  a  minimum  cathode  thickness 
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of  ~20  to  30  microns  ensures  the  applicability  of  Eq.  (37). 
A  similar  equation  is  also  applicable  for  the  anode.  For  the 
calculations  presented  here,  therefore,  Rcct  and  R^t  corre¬ 
spond  to  composite  electrodes. 

For  the  present  calculations,  the  charge  transfer  reactions 
are  assumed  to  occur  at  the  physically  distinct  cathode- 
electrolyte  and  anode-electrolyte  interfaces.  This,  of  course, 
is  not  strictly  correct  since  the  very  nature  of  composite 
electrodes  allows  for  the  occurrence  of  charge  transfer 
reactions  throughout  the  electrode.  In  this  simplified,  but 
analytical  approach,  however,  it  is  assumed  to  occur  at  the 
physically  distinct  electrolyte/electrode  interface.  This  is 
justified  as  long  as  X  is  much  less  than  tc. 

4.  Conclusions 

The  effects  of  interconnect  contact  spacing  and  contact 
area,  in  conjunction  with  transport  properties  and  design 
parameters  of  a  single-cell  on  overall  SOFC  stack  resistance 
were  investigated.  Area-specific  resistance  (ASR)  of  cell- 
interconnect  repeat  unit  was  used  as  the  figure  of  merit  for 
the  purposes  of  comparison.  Analytical  expressions  for  the 
ASR  for  interconnect  geometries  with  one-  and  two-dimen¬ 
sional  symmetries  were  derived  using  a  ladder  network 
approach.  The  advantage  of  analytical  expressions  over 
numerical  models  is  that  it  is  in  principle  possible  to 
explicitly  evaluate  the  role  of  various  parameters  on  the 
overall  stack  resistance  to  assist  in  the  optimization  of  stack 
performance.  The  analytical  results  then  may  guide  a 
detailed  numerical  modeling  effort.  While  the  one-dimen¬ 
sional  symmetry  is  expected  to  give  a  lower  resistance  than 
the  two-dimensional  symmetry,  the  present  analysis  gives 
closed  form  analytical  expressions,  which  allow  for  a  quan¬ 
titative  comparison  between  the  two.  For  an  interconnect 
design  with  a  one-dimensional  symmetry,  the  area-specific 
resistance  (ASR)  was  found  to  increase  linearly  as  a  function 
of  the  interconnect  spacing  for  widely  spaced  interconnect 
contacts,  and  the  size  of  the  interconnect  contact  area  had 
only  a  minor  effect.  The  area-specific  resistance  (ASR)  for 
an  interconnect  design  with  a  two-dimensional  symmetry, 
on  the  other  hand,  was  found  to  increase  parabolically  as  a 
function  of  interconnect  contact  spacing  for  widely  spaced 
interconnects  and  was  found  to  be  singular  (infinite)  at  zero 
interconnect  contact  area.  Thus,  interconnects  with  widely 
spaced  contacts  of  a  small  contact  area  should  be  avoided 
when  an  interconnect  design  when  a  two-dimensional 


symmetry  is  used.  Two  designs  of  single  cells  were  con¬ 
sidered  and  compared,  electrolyte-  and  anode- supported. 
Regardless  of  the  interconnect  design  symmetry  or  whether 
electrolyte-  or  anode- supported  cells  are  used  to  build  a 
SOFC  stack,  the  area-specific  resistance  (ASR)  begins  to 
rise  rapidly  due  to  sheet  resistance  within  the  electrodes  for 
interconnect  contact  spacing  greater  than  5  mm  for  the 
values  of  parameters  used  in  the  present  calculations.  The 
best  design  was  found  to  be  the  one  having  anode- supported 
cells  with  closely  spaced  interconnect  contacts  and  a  one¬ 
dimensional  symmetry. 

Acknowledgements 

This  work  was  supported  by  the  Gas  Research  Institute 
under  contract  #50942603276. 


References 

[1]  S.  de  Souza,  S.J.  Visco,  L.C.  De  Jonghe,  Reduced-temperature  solid 
oxide  fuel  cell  based  on  YSZ  thin-film  electrolyte,  J.  Electrochem. 
Soc.  141  (3)  (1997)  L35-L37. 

[2]  J.W.  Kim,  A.V.  Virkar,  K.Z.  Fung,  K.  Mehta,  S.C.  Singhal, 
Polarization  effects  in  intermediate  temperature  anode-supported 
solid  oxide  fuel  cell,  J.  Electrochem.  Soc.  146  (1)  (1999)  69-78. 

[3]  A.V.  Virkar,  J.  Chen,  C.W.  Tanner,  J.W.  Kim,  The  role  of  electrode 
microstructure  on  activation  and  concentration  polarizations  in  solid 
oxide  fuel  cells,  Solid  State  Ionics  131  (1-2)  (2000)  189-198. 

[4]  T.  Kenjo,  S.  Osawa,  K.  Fujikawa,  High-temperature  air  cathodes 
containing  ion  conductive  oxides,  J.  Electrochem.  Soc.  138  (2) 
(1991)  349-355. 

[5]  T.  Kenjo,  M.  Nishiya,  LaMn03  air  cathodes  containing  Zr02 
electrolyte  for  high-temperature  solid  oxide  fuel  cells,  Solid  State 
Ionics  57  (3-4)  (1992)  295-302. 

[6]  C.W.  Tanner,  K.Z.  Fung,  A.V.  Virkar,  The  effect  of  porous  composite 
electrode  structure  on  solid  oxide  fuel  cell  performance.  I. 
Theoretical  analysis,  J.  Electrochem.  Soc.  144  (1)  (1997)  21-30. 

[7]  E.F.  Sverdrup,  C.J.  Warde,  R.L.  Eback,  Design  of  high-temperature 
solid-electrolyte  fuel  cell  batteries  for  maximum  power  output  per 
unit  volume,  Energy  Conversion  13  (1973)  129. 

[8]  S.  Ahmed,  C.  McPheeters,  R.  Kumar,  Thermal-hydraulic  model  of  a 
monolithic  solid  oxide  fuel  cell,  J.  Electrochem.  Soc.  138  (9)  (1991) 
2712-2718. 

[9]  N.F.  Bessette,  W.J.  Wepfer,  J.  Winnick,  A  mathematical  model  of  a 
solid  oxide  fuel  cell,  J.  Electrochem.  Soc.  142  (11)  (1995)  3792- 
3800. 

[10]  C.R.  Wylie,  Advanced  Engineering  Mathematics,  McGraw-Hill, 
New  York,  1975,  p.  403. 

[11]  M.  Abramowitz,  I. A.  Stegun,  Handbook  of  Mathematical  Functions, 
Dover,  New  York,  1972,  p.  374. 


